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Abstract. We show that the introduction of a hyperbolic phase for Brans-Dicke 
(BD) field results in a flat vacuum cosmological solution of Hubble parameter H and 
fractional rate of change of BD scalar field, F which asymptotically approach constant 
values. At late stages, hyperbolic phase of BD field behaves like dark matter. 
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1. Introduction 

It has always been one of the most challenging and interesting problems of cosmology 
what the composition of the universe exactly is: what was it in the primordial time 
and what is it in today's universe? Where did the structure of the universe originally 
come from? After the development of the inflationary theory both observational 
and theoretical studies have been continuing on this subject. According to inflationary 
Universe models [2], inflation is capable of explaining not only the acceleration of 
the expansion rate but also flatness, homogeneity and isotropy of the universe. In 
addition, the discovery of the cosmic microwave background [3] indicates that our 
universe is nearly fiat and expands with a slow accelerating rate [1]-[10]. This slow 
rate acceleration of universe results from an adequate negative pressure of dark energy 
and recent observations indicate that dark energy behaves like Einstein's cosmological 
constant [11] which arises from the vacuum energy. The remaining energy density is 
composed of dark matter which can not be observed directly although its gravitational 
effects on visible matter validate its presence. In respect of recent WMAP data [12], our 
universe is composed of 72 % dark energy, 23 % dark matter, and 5 % ordinary (visible) 
matter. 

Up to now, the most popular candidate of the dark energy is the cosmological 
constant (vacuum energy) with the equation of state parameter u = —1. However, 
the observed vacuum energy density is at least 120 orders of magnitude smaller than 
predicted by particle physics. This is the so-called cosmological constant problem. In 
order to solve this problem, alternative models based on a dynamical cosmological 
constant A, with a negative equation of state have been constructed. These models 
include a scalar field with a slowly varying energy density. In quintessence models, the 
scalar field which is minimally coupled to gravity with an equation of state u > —1 
acts as dark energy and a potential energy dominating over kinetic energy leads to 
the accelerating expansion [13]. If the scalar field has a non-canonical kinetic energy 
then we have k-essence models [H]. On the other hand, phantom energy models with 
a negative kinetic energy assert an equation of state parameter u < —1 [13]. Besides, 
string-theory inspired quintom models have also been analyzed. A model which includes 
the combination of two-scalar fields have been considered [15] , one corresponding for the 
early time quintessence dominance, uj > —1 and the other one corresponding for the 
late time dominance, u; < — 1. In addition, another string inspired quintom model 
where tachyon is non-minimally coupled to gravity obtained the conditions required 
for UJ crosses over —1 [16]. Modified gravity models in the framework of scalar-tensor 
theories have also been analyzed to explain the acceleration of the universe [IT]. A 
special case of these type of models is the Brans-Dicke- Jordan- Thirry [l8]-[20] theory 
where the curvature scalar occurs only linearly in the lagrangian density. Whether 
the quintessence field can be identified with the Brans-Dicke- Jordan-Thirry field is an 
interesting question [2T]-[25]. In addition to explaining dark matter, BD theory may 
have other advantages. In particular it has been remarked that BD theory can be 
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imbedded in electroweak theory [22] and it can explain the cosmic coincidence problem 
[30] . There exist a number of studies on accelerated models in BD theory [3T|- [37] . For 
example, Sen et al [38| have found the potential relevant to power law expansion in BD 
cosmology. In addition, in a work of Setare [39|, the lower bound of cja was found —0.9 
using a holographic dark energy model in the framework of non-flat BD cosmology. 

In standard cosmology the rate of expansion of the universe strongly depends on 
the equation of state of the matter-energy that fills it. One immediate question which 
arises is that whether there is any consistent modification of Einstein's equations such 
that the expansion of the universe is independent of its content. In the previous works 
of Arik, Calik and Sheftel [23|-[25|, it is shown that BD scalar tensor theory of gravity 
with the standard mass term potential (l/2)m^0^ is capable of explaining the rapid 
primordial and slow late-time inflation and a linearized non-vacuum late time solution 
well accounts for the contribution of dark energy to the Friedmann Equation, however, 
it does not account for the contribution of dark matter. In this regard, we particularly 
focus on the model consists of a modified Brans-Dicke- Jordan-Thirry [IE]- [20] model 
where both the signs of the kinetic term {(fxp* = — (p^) in its 02 part and potential 
term bring a minus sign. The models with this sign convention in Lagrangian have been 
termed as quintom models [15], [16], a word induced from quintessence and phantom. 
We add an imaginary part to the BD field such as = 0i + i02, and search for a 
contribution to dark matter in the presence of the imaginary part of (p field, 02- 

2. Field Equations 

In this work, we will show that both the dark matter contribution ^dm and dark 
energy contribution to Friedmann Equation can be explained solely by BD theory 
of gravity provided that BD scalar field is modified suitably. The most straightforward 
modification is choosing the BD field as a complex field defined by 



where 0r is real scalar field amphtude. Such complex BD field can also be represented 
as in matrix form 



where (T2 is a Pauli spin matrix. 

However, we will take the phase of to be hyperbolic by replacing the term z/3 = \l/ 
in ([1]) such that becomes 



'1 + 202 = 0Re' 



(1) 




(2) 




(3) 



and its conjugate matrix becomes 




(4) 
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where 

01 = cosh ^ (5) 

4>2 = 4>R sinh (6) 

where \1/ is reaL With this modification, we note here that gains a "Quintom" 
character since its kinetic contribution {(fxp* = (pl — (pl) brings a minus sign. But 
nevertheless, in this paper, we will show that a cosmological vacuum solution with flat 
space-like section is capable of explaining how the Hubble parameter H evolves with the 
scale size of the universe a{t) and how the solution of fractional rate of change of BD 
scalar field, F contributes to the evolution of H in the late era in which the universe is 
expanding at a slow rate. In the context of BD theory [20] with self interacting potential 
and matter field, the action in the canonical form for real BD scalar field in our notation 
is given by 



oUJ Z Z 



(7) 



however, since we have modified the scalar BD field as in (jll), we also modify the 
action above as 

'-1 , 1 1 



S" = ^ tr y d^'x^ 



(8) 



where / is the unit matrix. In particular we may expect that is spatially uniform, 
but varies slowly with time. The signs of the non-minimal coupling term and the 

kinetic energy term are properly adopted to (H ) metric signature. In units where 

c = = 1, we define Planck-length, Lp, in such a way that L\(f)\ = uj/2it where (pn is 
the present value in (I5f6il . Hence the dimension of the scalar field is chosen to be L~^, so 
that Geff has a dimension Lp since nonminimal coupling term R where R is the Ricci 



Tgjj- nab a umieusiou ijp binue uouimmiiim coupimg leim ipj^. 

scalar, replaces with the Einstein-Hilbert term -^R in such a way that G~jj = ^(pR 
where Ge/f is the effective gravitational constant as long as the dynamical scalar field 
varies slowly with time. To be in accordance with the weak equivalence principle, the 
matter part of the Lagrangian, Lm, is decoupled from such that we have considered 
the energy-momentum tensor = diag (p, —p, —p, —p) just with classical perfect fluid 
where p is the energy density, p is the pressure. The gravitational field equations derived 
from the variation of the action ([8]) with respect to Robertson- Walker metric is 
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-^ + '^ - + - + - =0 1^ 

a J (p 2uj \a a^y 

where k is the curvature parameter with k = —1, 0, 1 corresponding to open, flat, 
closed universes respectively and a (t) is the scale factor of the universe (dot denotes 
^). Since in the standard theory of gravitation, the total energy density p is assumed 
to be composed of p = + Pm where pA is the energy density of the universe due to 
the cosmological constant which in modern terminology is called as ^^dark energy", the 
right hand sides of ([91 [H]) are adopted to the matter energy density term pM instead of 
p and pm instead of p where M denotes everything except the field. The main reason 
behind doing such an organization is that whether if the terms on the left-hand side of 
([9]) can accommodate a contribution to due to what is called dark matter. In addition, 
the right hand side of the equation (fTTj) is set to be zero according to the assumption 
imposed on the matter Lagrangian Lm being independent of the scalar field 0. For the 
vacuum {p = p = 0) and fiat space like {k = 0) section solutions, we start with defining 
the fractional rate of change of as 

where 

0101 -0202 0R 0201 -0102 ^j, , 

n - 02 (i>R n - n 

and the Hubble parameter as H (a) = a/a, hence, we rewrite the left hand-side of the 
field equations (P fTTj) in terms of H{a), Fi{a), F2{a) and their derivatives with respect 
to the scale size of an universe a, as 

3H^-2uF^+2uF^+6FiH-2um'^ = (14) 

3H^ + {2u + 4) - 2cuF| + AFiH + 2aHF[ + 2aHH' - 2ujm^ = (15) 

- 6//^ + 2uF^+2ujF^ + 6uFiH+2uaHF[ - 3aHH' + 2ujm^ = (16) 

{AujFi+QujH) F2 + 2ujaHF2 = (17) 

where prime denotes Since solving these coupled equations analytically is hard 
enough, we have put forward following perturbation solution as 



F2 = f2o. + f2l(^)%F22(^)'° (20) 

where Hoo, Hi, H2, -Fioo, -^11, -^12, -^200, -P21, -F22 are perturbation constants and a is an 
exponential factor to be determined. With the transformation 
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3H^ - 2uF^ + 2ujF^ + QHFi - 2ujm^ = (22) 

3H^ + {2uj + 4) - 2cjF| + AHFi - 2auH ( ^ ) - 2auH ( ^ ) - 2ujm'^ = (23) 



du J \ du J 

- + 2ujF( + 2ujF'i + QuHFi - 2ujauH ( ^ ) + 3auH ( ^ ) + 2ujm^ = (24) 



du J \du J 

- 2uauH (^^^ + {^ujFi+QujH) F2 = (25) 

and f|T8ll20l) becomes 

H = H^ + Hiu + H2U^ (26) 

Fi = Flo, + i^iiM + i"i2M' (27) 

F2 = F200 + F21U + F22U'. (28) 

Substituting fl26ti28p into fl22ti25p and keeping only the zeroth, first and second order 
terms of u and neglecting fiiglier order terms of m, we get tlie following equations to be 
solved. In the zeroth order of u; 

3Hl-2ujFl^+2u:Fl^+QF^^H^-2ujm'' = (29) 

ml, + {2uj + 4) Fl^ - 2uoFl^ + 4Fiooi^oo - 2um^ = (30) 

- 6Hl + 2uoFl^+2ujFl^ + QujFi^H^+2tum'' = (31) 
[4cuFioo+6cuifoo]F2oo = (32) 

in the first order of u; 

(6F100 + QHoo)H, + (6//oo-4cuFioo)Fn + 4ujF2iF2oo = (33) 

{{6-2a)H^ + AFi^)H^ + {{A-2a)H^ + {Au + 8) Fioo)Fn - AU0F21F200 = (34) 

((3a - 12)H^ + QuoFi^)Hi + ((6cu-2cja)/7oo + AujFi^)Fii + AU0F21F20C = (35) 

[(-2cja + 6cj)i7oo + AujFi^]F2i + ALuFnF2oc + Qu)H^F2oo = (36) 
in the second order of m; 

ml + 6i/oo^2 - 2cuF2, - 4cuFiooFi2 (37) 
+4CUF200F22 + 2cuF|i + 6Fi„oi^2 + 6Fn//i + 6F12//00 = 

(3 - 2a)Hl+ (2cu+4) F^'i + (4 - 2a)FuHi + {AF,^ + 6i^oo - AaH^)H2 (38) 
+ [4iyoo + (4cu+8) Fioo - AaH^]Fi2 - 2cuF|i - 4c^F2ooF22 = 



(3a - Q)Hl+2uFl^ + (-2cja + Quj)FiiHi + (-12i7oo + QaH^ + QuFi^)H2 (39) 
+ (-4cua^roo + 4cuFioo + 6ujH^)Fi2 + 2cjF|i + AUJF200F22 = 



(4a;Fioo - 4a;aiyoo + QujH^)F22 + {-2ujaHi + AuFn + 6ljHi)F2i (40) 
+4cx;F2ooFi2 + 6UJF200H2 = 0. 
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Solving the equation set ( 12911321) and (13311361) provide respectively, 

F2^ = H^= + F,^ = J^ (41) 

v/(6cj2 + ij^ + 12) 2c^ + 2 

1 3 

a = 3 H Fii = — Hi F21 = free — parameter (42) 

CO + 1 2 

and afterwards, substituting (HT| H2|) into the equation set (13711401) yields the following 
equation set to be solved for H2, F12, F21, F22 as 

(12cj + 18)H^H2 + {8uj + 12)H^Fi2 + 4cj(cj + l)F|i = {9uj^ + 2luj + l2)Hl (43) 
- {l2ujlQ)H^H2 - (12cj + lQ)H^Fi2 - 4cu(cj + 1)^2^ = -(9cj2 + 27tu + 2Qi)Hl (44) 
(ISc^ + 2A)H^H2 - (12cj^ + 16c<;)i7ooi^i2 + 4cu(cu + 1)F|i = -(Qcj^ + 21cu + 12)/72 (45) 

^'-F2i. (46) 



22 



00 



To proceed one step further, we write the standard Friedmann equation: 

|)^«..«„(^)' (47, 

and we fit all theory parameters to the observational density parameters; 

= (48) 

= (49) 



where 



Hl = Hl + 2H^{Hi + H2) + Hi (50) 



With these relations above and the constraint + Qm = 1, where Qm = ^vm + 
^DM; we can express theoretical parameters Hi in terms of the observational density 
parameters Qa, and H^o as 

Using recent observational results [12j on density parameters Qdm — 0.28, Qa ^ 0.72 
and fivM = (since the universe we study in this theory is vacuum) together with (15Ti) 
we determine; 

Hi = ^H^^0.19H^. (52) 
Similarly, when we solve the equations (1431146 P : the solutions are; 

H2H0C, = (iSu^H^ - Su^Fl - Au'^F^, - 12H^ - uoH^ - AuF^, + 9u^Hl;) (53) 
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As — > 00: 



'12 



3 i^/^oo 4 i^ooy 

^ -^li + TT^i^n (56) 



24if^ 2AH, 



00 



i^22 = -^i"21 (57) 
-noo 

At this point, we emphasize that H2 must be zero in order to make sense with 
recent observational data on density parameters of the universe and to find the exact 
value for F21. Therefore, we insert H2 = and we get 

F21 ~ (58) 

F12 - OmH^ (59) 

F22 ^ -OMHoo. (60) 

Hence, with these perturbation constants fHT| H2| [52l [58]|60|) found from theory we 
can express ( fT8ll20l) : 

H = H^ + 0.19H^(^^y (61) 

F, = - 0.28H^ (^Y + OmH^ (62) 

2^1; + 2 \ a / \ a / 

r2 = U.ZOiloo I 

where 



F2 = 0.28i/oo(-) -0.06ifoo(^) (63) 



~ 0.84i/o (64) 
if (16T!) is satisfied for = i^o, and i^o is the present value of the Hubble parameter [T2] . 



4. Conclusion 

In this paper, we have analyzed the dark matter Qdm and dark energy contribution 
to Friedmann Equation solely by modified BD theory of gravitation with no other 
input. As far as we know, the scalar field was always examined individually, however, 
we brought forward a new idea such that it can have different components and each of 
these components can account for different energy densities. 

Actually, the starting point of our motivation originates from this point in the sense 
that when BD theory of gravitation with solely scalar field (p is substituted into role as 
discussed in our previous work [25j, we have shown that WMAP+SnIa data [12], [27]-[29] 
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favor this model instead of the standard Einstein cosmological model with cosmological 
constant (LCDM model [2H]-[29]) under the condition that the new density parameter 
fl^ induced in Friedmann equation in standard cosmology to be < and H2 seen in 
the equation (l26l) to he H2 < instead of H2 = 0. (for further information, see [25]). In 
other words, at this stage, we have realized that the more we force H2 to be less than zero 
in the real phase of the model with individual scalar field 0, it fits WMAP+SnIa data 
much more confidentially than LCDM model |28]-[29]. To do this, in the first attempt, 
we have used a complex scalar field = 0i + i(f)2 so that WMAP+SnIa [I2], [27]-[29] 
data will favor the model with modified BD field in its complex phase. Although this 
approach has brought brand new considerations and aspects to Friedmann Equation in 
the concept of dark matter and dark energy, a more suitable solution was found with 
the modification of scalar field by using a hyperbolic phase i(3 = . Having solved the 
field equations including the hyperbohc phase, we achieved the field equations ( |T4lfT7l) 
of modified BD scalar tensor theory namely equations of "Quintom" model. 

To solve these above field equations, we put up the argument of perturbative 
solutions with the constant terms H^o, Hi, H2, Fi^oFu, F12, -F2005 -^21 and F22. All of these 
constants have made it possible to originate new predictions on dark matter and dark 
energy contribution of BD theory. 

To begin with, the most significant evidence for the idea that this modification 
needs real attention is the solution of a. It can easily be seen that, when u 00 (where 
BD approaches Einstein theory), a —> 3, as it appears in the Friedman Equation in the 
form Qm (^)^- Similarly, the term H^o which has no scale factor term, just like the 
energy density term due to the cosmological constant Q/^ in the Friedman Equation, 
was found purely from theory; 

v/(6cu2 + + 12) 

From the equation (16T1) in the equation set fl6m63p , we see that the second term is found 
to be smaller than the first one and the third term is found to be smaller than the 
second one. Namely, the dominating term is the first one which can be interpreted as 
the contribution to dark energy. On the other hand, the second term can be considered 
as the contribution to dark matter. However, the situation is different for Fi and F2 as 
it is seen in the equations (162| [63]) . As it was mentioned before, in the absence of F2 
term, where F = Fi, the theory was able to explain the contribution to dark energy but 
not to dark matter. Our aim was to find a contribution to dark matter in the presence 
of F2, with the component F21 since it is coupled with {^)^- Namely; it is agreeable to 
predict that while Fi^o which is not coupled with a scale factor term is contributing to 
dark energy, F21 is contributing to dark matter. Hence, the introduction of a hyperbolic 
phase for BD field results in a flat vacuum cosmological solution of Hubble parameter 
H and fractional rate of change of BD scalar field, F which asymptotically approach 
constant values. At late stages, hyperbolic phase of BD field behaves like dark matter. 
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